Scalable Parallel SSOR Preconditioning for Lattice Computations in Gauce Theories
نویسندگان
چکیده
We discuss a parallelization scheme for SSOR precondition-ing of Krylov subspace solvers as applied in lattice gauge theory computations. Our preconditioner is based on a locally lexicographic ordering of the lattice points leading to a parallelism adapted to the parallel system's size. By exploitation of thèEisenstat-trick' within the bi-conjugate gradient stabilized iterative solver, we achieve a gain factor of about 2 in the number of iterations compared to conventional state-of-the-art odd-even preconditioning. We describe the implementation of the scheme on the APE100/Quadrics SIMD parallel computer in the realistic setting of a large scale lattice quantum chromodynamics simulation.
منابع مشابه
A Scalable Parallel SSOR Preconditioner for Efficient Lattice Computations in Gauge Theories
solver We extend the parallel SSOR procedure for the eecient preconditioning of modern Krylov subspace solvers 1], recently introduced in 2] towards higher order, quantum-improved discretization schemes 3] for lattice quantum chromodynamics 4]. 1. Solving linear systems in LGT. Lattice gauge theory (LGT) deals with the controlled numerical evaluation of gauge theories like quantum chromodynamic...
متن کاملParallel SSOR preconditioning for lattice QCD
The locally-lexicographic symmetric successive overrelaxation algorithm (ll-SSOR) is the most eeective parallel preconditioner known for iterative solvers used in lattice gauge theory. After reviewing the basic properties of ll-SSOR, the focus of this contribution is put on its parallel aspects: the administrative overhead of the parallel implementation of ll-SSOR, which is due to many conditio...
متن کاملPreconditioning of Improved and “Perfect” Fermion Actions
We construct a locally-lexicographic SSOR preconditioner to accelerate the parallel iterative solution of linear systems of equations for two improved discretizations of lattice fermions: (i) the Sheikholeslami-Wohlert scheme where a non-constant blockdiagonal term is added to the Wilson fermion matrix and (ii) renormalization group improved actions which incorporate couplings beyond nearest ne...
متن کاملA Preconditioner for Improved Fermion Actions
SSOR preconditioning of fermion matrix inversions which is parallelized using a locally-lexicographic lattice sub-division has been shown to be very efficient for standard Wilson fermions. We demonstrate here the power of this method for the Sheikholeslami-Wohlert improved fermion action and for a renormalization group improved action incorporating couplings of the lattice fermion fields up to ...
متن کاملParallel Implementation of Preconditioned Iterative Schemes by Means of Overlapping Decomposition
We consider iterative solvers for large linear systems, and develop a theoretical analysis for parallel implementation techniques which resort to overlapping decompositions, i.e. in which some unknowns are replicated on two or more processors to facilitate the parallelization. We investigate the conditions under which a given sequential preconditioner can be implemented by means of purely local...
متن کامل